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Abstract 

We present a prescription to compute the time-domain gravita¬ 
tional wave (GW) polarization states associated with spinning compact 
binaries inspiraling along quasi-circular orbits. We invoke the orbital 
angular momentum L rather than its Newtonian counterpart Tn to 
describe the orbits and the two spin vectors are freely specified in the 
source frame associated with the initial direction of the total angular 
momentum. We discuss the various implications of our approach. 


1 Introduction 

Gravitational waves (GWs) from coalescing compact binaries containing at 
least one spinning component are expected to be routinely detected by 
the second-generation laser interferometric detectors like advanced LIGO 
(aLIGO), Virgo and KAGRA [T]. The detection of GWs from such bina¬ 
ries and subsequent source characterization crucially depend on accurately 
modeling temporally evolving GW polarization states, h^(t) and hx (t), from 
such binaries during their inspiral phase [2]. At present, and hx{t) as¬ 

sociated with non-spinning compact binaries inspiraling along quasi-circular 
orbits have GW phase evolution accurate to 3.5PN order and amplitude 
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corrections that are 3PN accurate [3]. Recall that the 3.5PN and 3PN or¬ 
ders correspond to corrections that are accurate to relative orders {v/cY 
and (u/c)® beyond the ‘Newtonian’ estimates, where v and c are the orbital 
and light speeds, respectively. In the case of inspiraling compact binaries 
containing Kerr black Holes (BHs), it is desirable to employ temporally 
evolving /i+(t) and hx{t) that incorporate the spin effects very accurately 
and the dominant spin effect arises due to the general relativistic spin-orbit 
coupling that appear at the relative 1.5PN order for maximally spinning 
Kerr BHs BM- Following Ref. [S], we define the spin of a compact object 
as S = GrrY^x^lcj where mco,x and s are its mass, Kerr parameter and 
a unit vector along S, respectively and for a maximally spinning Kerr BH 
X = 1. We note that it is customary to employ the Newtonian orbital an¬ 
gular momentum Rn = fir x v, where fi,r and v are the reduced mass, 
orbital separation and velocity, respectively, to specify these quasi-circular 
orbits [SI [6]. 

In what follows we provide a prescription to generate the time-domain 
amplitude corrected /i+(t) and hx{t) for spinning compact binaries inspi¬ 
raling along quasi-circular orbits, described by their orbital angular mo¬ 
menta. Note that the amplitude corrected /i+(t) and hx{t) refer to GW 
polarization states that are PN-accurate both in its amplitude and phase. 
We describe our approach to compute the fully I.5PN accurate amplitude 
corrected expression for hx (t) while invoking L, the PN-accurate orbital an¬ 
gular momentum, to characterize the binary orbits. We symbolically obtain 
the additional I.5PN order amplitude corrections to hx{t) in comparison 
with equations (A3) in Ref. |6] that employ Ln to describe the binary or¬ 
bits. We also discuss certain implications of our approach while considering 
spin effects due to the leading order general relativistic spin-orbit coupling. 
Our attempt to perform GW phasing with the help of L is motivated by a 
number of observations (we term accurate modeling of temporally evolving 
GW polarization states as ‘GW phasing’). First being Ref. [7] that provided 
a prescription to implement GW phasing for spinning compact binaries in 
inspiralling eccentric orbits in an accurate and efficient way. We are fur¬ 
ther influenced by the fact that it is customary to use precessional equation 
appropriate for L to evolve Rn while incorporating the effects due to the 
dominant order spin-orbit coupling [SIEIE]. Finally, we note that a seminal 
paper that explored the inspiral dynamics of spinning compact binaries and 
the influences of precessional dynamics on /i+^x(i) employed L to describe 
their binary orbits [9]. Additionally, we specify the two spins in an inertial 
frame associated with the initial direction of the total angular momentum 
Jq while it is customary to invoke a non-inertial LN-based orbital triad to 
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specify the two spins in the literature. 


2 GW phasing for spinning binaries characterized 
by L 

We invoke, as noted earlier, the PN-accurate orbital angular momentum L 
to describe binary orbits and we specify at the initial epoch both the orbital 
and spin angular momentum vectors in an inertial frame associated with 
the initial direction of total angular momentum (see figured]). We begin by 
presenting an expression for the cross polarization state having Newtonian 
(quadrupolar) order amplitude and the relevant expression reads [TU] 


^xlQ(f) = 2 


Gn 

c2 R' ^ 


(1 


cos i) Sq sin L sin(Q; — 2<1>) 


— (1 + cos i) Sg sin i sin(a + 2<h) 

— -(1 + 2 cos i + cos^ i) Cg sin(2Q; + 2<h) 

— - (1 — 2 cos L + cos^ l) Cg sin(2a — 2<f>) > , 


( 1 ) 


where R',Sg and Cg stand for the radial distance to the binary, sin0 and 
cos 9, respectively. The dynamical angular variable $ measures the orbital 
phase from the direction of ascending node in a plane perpendicular to k 
while its derivative is required to define jc? = (Gm To obtain 

temporally evolving hxlqit) associated with spinning compact binaries in- 
spiraling along quasi-circular orbits, we pursue the following steps. First, 
we specify how the Eulerian angles a and l that specify the orientation of 
L vary under the conservative orbital dynamics. This requires us to employ 
three differential equations that describe the precessional dynamics of the 
orbital and spin angular momenta and these differential equations contain 
an orbital like frequency uj [5|. The conservative evolution for $ and d* are 
governed by the following differential equation 

= w — cos bOi, (2) 

where the orbital like frequency ui is defined by the relation ui = vjr 
that connects uj to the orbital separation and velocity. This PN-accurate 
equation arises from the vectorial expression for the orbital velocity in 
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Figure 1: The inertial frame and the Cartesian coordinate system where 
the z axis points along Jq, the direction of total angular momentum at the 
initial epoch. We display the angles that characterize the orbital and spin 
angular momentum vectors, denoted by k, si and S 2 , while the line of sight 
vector N is in the x — z plane. The dashed lines depict the projections of 
various vectors on the x — y plane. 


a CO- moving triad [5]. It is convenient to introduce a dimensionless pa¬ 
rameter X = {Gmuiorh/c^)‘^^^ such that the above equation reads ‘h = 
(x^/^/(Gm/c^)) — cos id. Thereafter, we impose the effects of gravitational 
radiation reaction on these differential equations by specifying how x, ap¬ 
pearing in these differential equations, vary during the binary inspiral. 

It turns out that we require to solve simultaneously the precessional 
equations for fc,si and S 2 , the unit vectors along L,Si and S 2 , to spec¬ 
ify how OL and i vary under the conservative dynamics as the precessional 
dynamics of L^S\ and S 2 are intertwined. In figure [H we display various 
angles that specify the orientations of these three unit vectors in the above 
inertial frame. The Cartesian components of these vectors are 

k = (sin i cos a , sin i sin a , cos C ), (3a) 

51 = (sin 01 cos (/)!, sin 01 sin (/)i, cos 0i) , (3b) 

5 2 = (sin02 cos(/> 2 , sin02 sin(/) 2 , cos 02) . (3c) 

The precessional equations for the above three unit vectors, extractable from 
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Ref. [HIS], read 



(4b) 


(4a) 


(4c) 


where q = mi/m 2 and the symmetric mass ratio rj = n/m is required 
to define the quantities (5i and <52 ( hi ,2 = rj/2 + 3(1 =F ~ 4r/)/4). We 
are now in a position to incorporate the effect of gravitational radiation 
damping and this is achieved by specifying the secular variations in x in 
the above differential equations. The PN-accurate differential equation for 
X arises from the usual energy balance arguments and may be obtained from 
equation (3.21) in Ref. [6]. 

We obtain numerically the time-domain GW polarization states for our 
binaries by simultaneously solving the differential equations for the Carte¬ 
sian components of k,si and S 2 along with the PN-accurate equations for 
$ and dx/dt. This implies that we numerically solve a set of 11 differential 
equations and the resulting variations in a, l, and $ are implemented in 
equation ([T]) to obtain the temporally varying hxlqit). A close inspection 
reveals that we require to specify four angles that provide the orientations 
of the two spins from j'q at the initial epoch to obtain hxlqit). It should 
be noted that the two angles that specify the initial orientation of k from 
Jq are not independent variables. This is because one may point the total 
angular momentum along the 2 ;-axis at the initial epoch without the loss 
of any generality. This allowed us to equate the x and y components of 
J = L + Si + S 2 to zero at that epoch and therefore to obtain the initial 
estimates for l and a. In our numerical integrations, the bounding values 
for x are given by xq = 2.9 x 10“^ {m' and Xf = 1/6, where m' is the 

total mass of the binary in solar units and let cjq = IOtt Hz as customary 
for aLIGO along with <h(xo) = 0. Finally, we note that the values of a and 
i at every step of our numerical runs are obtained from the Cartesian com¬ 
ponents of k by the relations: a = cos“^(A:x/y^fe|4-^) and l = cos“^(fcz). 
Therefore, we require eight independent parameters to specify the inspiral 
dynamics of compact binaries with spinning components and these eight pa¬ 
rameters are the four basic (constant) parameters namely (m-i, m 2 , Xi; X 2 ) 
along with the above four angular dynamical variables that specify the two 
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spin vectors in the inertial frame, namely and ( 02 ,</' 2 ) as displayed 

in figure [T] 

In what follows we compare our approach with what is detailed in Ref. [6] 
that provided a way to obtain the time-domain ready-to-use GW polariza¬ 
tion states for inspiraling compact binaries while incorporating all 1.5PN 
order spin effects both in amplitude and GW phase evolutions. Their ap¬ 
proach differs from ours in two aspects. First, Ref. [6] employed I = Ln/\ 
to characterize the orbital plane and following the earlier papers employed a 
precessional equation for I that is identical to our equation (I4ap for k while 
replacing fc by Z. It is not very difficult to show that this is equivalent of 
using an orbital averaged expression for I [5]. However, invoking an orbital 
averaged precessional equation for I leads to an undesirable feature that the 
coefficient of I in the expression for h in the (n, X = I x n,l) frame will not, 
in general, vanish. It is fairly straightforward to compute the time derivative 
of n and express it in the (n. A, Z) frame as 


dn 

dt 


fd^' 

— -h COS i 

\ dt 


/ 


da'\ 
dt ) 


A + 


f - I ^ I 

— sm $ — sm i cos —— Z, 
\dt dt) 


(5) 


where d, a' and are the usual 3 Eulerian angles required to dehne n, 
A and Z in the jg-based inertial frame. The fact that Z = n x n/|n x n| 
clearly demands that the coefficient of Z in the above equation should be 
zero as also noted in Ref. [Uj. However, if one employs equation (Hal) to 
describe the precessional dynamics of Z then it is possible to show with some 
straightforward algebra that 


, di , , da 

sm cp —-cos sm i —— = — —— 

dt dt G m 


hqxi («i ■ >d) + ^X2 {S2 ■ A) 


( 6 ) 


as noted in Ref. [7]. It is not very difficult to conclude that the right hand 
side of above expression, in general, is not zero. 

It turns out that h having components along Z leads to certain anomalous 
terms that contribute to the evolution at the 3PN order and this is within 
the consideration of higher order spin effects available in the literature. Note 
that the higher order spin effects are known to 3.5PN order while dealing 
with the spin-orbit interactions m- We demonstrate our observation by 
noting that the dehnitions v = rh and v = roj imply that uP = n fi. With 
the help of our Eq. ([5|) the expression for oj reads 
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(jj = 


1 


— -h COS i 

\ dt 


/ 


da''\ 
dt ) 


+ 


2 $' 


dd ,, , 

—— sin ® — sin L cos ® —— 
dt dt 


(7) 


where stands for d^'/dt and the higher order terms that are cubic in the 
time derivatives of d and a' are neglected. Invoking the fact that at the 
Newtonian order is given by /{Gm/(?) and noting that the expression 
for dd/dt and da'/dt arise from the 1.5PN order differential equation for k, 
we get 

(l + A' + , (8) 

Gm \ J 

such that A' and B' are given by 

dd = (si • A) + ^X 2 (S 2 • A)] cos^>' 

+ [(5i q xi (si • n) + ^ X 2 (^2 • n)] sin ^>'| , (9a) 

B' =-^ |(5igxi (si • A) + ^X2 (s2 • A)| . (9b) 

It should be evident that the B' terms appear at the 3PN order while A' 
terms enter <I> expression at the 1.5PN order. It is not very difficult to in¬ 
fer that these B' terms arise due to the non-vanishing I component in the 
expression for n and therefore are unphysical in nature. Therefore, these 
anomalous terms contribute to the <&' evolution at the third post-Newtonian 
order and this is within the consideration of higher order spin effects cur¬ 
rently available in the literature. It should be noted that these unphysical 
terms play no role in investigations in Ref. [6] that probed the leading order 
spin effects appearing at the 1.5PN order in the phase evolution. 

Another consequence of invoking k to specify the binary orbit is the 
appearance of certain new 1.5PN order contributions to the amplitudes of 
h+ and hx in addition to what is provided by equations (Al), (A2) and (A3) 
in Ref. [6]. We recall that equations (Al), (A2) and (A3) in Ref. [6] provide 
the fully 1.5PN accurate expressions for /i_|_ and hx while invoking Rn to 
describe the binary orbits. The additional amplitude corrections to the GW 
polarization states arise mainly due to the fact that the component of v along 
k is of 1.5PN order. To demonstrate this point, we express r and v = dr/dt 
in the inertial frame {x, y, z) associated with jig with the help of the three 
usual Eulerian angles <h, a and t as displayed in figure [H The relevant 
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expression for r may be written as r = rn where n = (—sina cos<^ — 
cos L cos a sin $)» + (cos a cos — cos i sin a sin ^)y + sin i sin ^z. To show 
that V can have non-vanishing 1.5PN order terms along k, we compute dr/dt 
in the co-moving frame defined by the triad (n, ^ = k x n,k) and this is 
easily achieved with the help of three rotations involving the three Eulerian 
angles appearing in the expression for r [12]. The resulting expression for v 
reads 


/ d^ da \ f di ^da\ 

\dt dt \dt dt J 


( 10 ) 


It is not very difficult to verify that v ■ k ^ 0 while invoking equation (|4al) 
for k to evaluate the k component of v. Moreover, the coefficient of k in the 
above expression for v is at the 1.5PN order. The k component of v enters 
the expressions for /i+ and hx through the dot products {p ■ v) {q ■ v) and 
{p ■ v)^ — (q ■ that are required to compute the PN-accurate expressions 
for /i+ and hx (the vectors p and q are two unit vectors that lay in a plane 
perpendicular to IV). The resulting 1.5PN order amplitude corrections to 
hx symbolically read 


hx 


1.5PN 


G fj, Gm 
R' y/x 


^ Ki(i, a, 4>, 6») -k ^ a, 4>, 9) 


( 11 ) 


where the explicit expressions for ki and K 2 are available in Ref. |10] . There¬ 
fore, the fully 1.5PN order amplitude corrected expression for hx associated 
with the spinning compact binaries in quasi-circular orbits, described by L, 
is provided by equations (A3) along with the above equation (fTTI) . This 
statement also requires that the angular variables i and a that appear in 
equations (A3) of Ref. [6] represent k rather than Z. Let us emphasize again 
that equations (A2) and (A3) in Ref. |6] indeed provide the fully 1.5PN 
accurate amplitude corrected and hx for spinning compact binaries in 
circular orbits described by Ln- 

Another aspect where we differ from Refs. lai is the way we specify 
the two spin and k vectors to perform GW phasing. In the literature, it 
is common to freely specify spin vectors in an orthonormal triad dehned 
by using I [SIISIE]. In contrast, we freely specify the two spin vectors at 
the initial epoch in the inertial source frame associated with j'q. This choice 
allowed us to specify the x and y components of k uniquely in terms of si, S 2 
and other intrinsic binary parameters at the initial epoch by demanding that 
the X and y components of J = L + Si + S 2 should be zero at that epoch. 


8 





Therefore, we extract very easily the initial values of a and i from the 
following expressions for the initial x and y components of k 


/cx,i = sin i cos a = -— {Xi xi sin cos (j)i + X 2 sin 62 cos (/> 2 } , 

C 

G m? 

ky^i = sin L sin a = --— {Xf xi sin sin + X| X 2 sin 02 sin ( 1 ) 2 } , 

C 

where 6 [, c/)}, 62 , <?i >2 the values of , <^i, ^2 j 4'2 at the initial epoch and Li 
denotes the PN-accurate expression for |L| at initial orbital frequency. 

For isolated unequal mass spinning compact binaries with g > 3 it should 
be advantageous to specify their spins in the inertial frame associated with 
Jq. This is because 9i{xq) that provides the dominant spin orientation from 
Jq at xq is expected to lie in a smaller range for such binaries spiraling into 
xq due to the emission of GWs. A recent study reveals that the dominant 
BH spin orientation from jig at xq is more likely to be < 90° for angular mo¬ 
mentum dominated unequal mass binaries (|L|(xo) > ^i) while 9i{xo) < 45° 
for unequal mass binaries having Si > \L\ at xq m- These inferences re¬ 
quire that these binaries inherit spin-orbit misalignments < 160° from var¬ 
ious astrophysical processes responsible for the formation of such unequal 
mass spinning binaries. The physical explanation for these conclusions is 
the observed alignment of Si towards jg due to the action of gravitational 
radiation reaction, detailed in Ref. [9]. However, it will be difficult to pro¬ 
vide similar bounds for the dominant spin orientation if the spins are freely 
specihed in a non-inertial orbital triad associated with L at xq- In this case 
k ■ Si provides the dominant spin orientation at xq and it is not difficult to 
show that k ■ si remains fairly constant as these unequal mass binaries with 
q > 3 spiral into xq from initial orbital separations < 1000 Gm/c^ due to 
the emission of GWs. 
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